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Abstract. We present an extensive study of magnetic fields in a system of merging galaxies. We obtained for 
NGC 4038/39 (the Antennae) radio total intensity and polarization maps at 8.44 GHz, 4.86 GHz and 1.49 GHz 
using the VLA in the C and D configurations. The galaxy pair possesses bright, extended radio emission filling 
the body of the whole system, with no dominant nuclear sources. The radio thermal fraction of NGC 4038/39 was 
found to be about 50% at 10.45 GHz, higher than in normal spirals. Most of the thermal emission is associated with 
star-forming regions, but only a part of these are weakly visible in the optical domain because of strong obscuration. 
The mean total magnetic fields in both galaxies are about two times stronger (~ 20 ^G) than in normal spirals. 
However, the degree of field regularity is rather low, implying tangling of the regular component in regions with 
interaction-enhanced star formation. Our data combined with those in Hi, Ha, X-rays and in far infrared allow 
us to study local interrelations between different gas phases and magnetic fields. We distinguish several radio- 
emitting regions with different physical properties and at various evolutionarystages: the rudimentary magnetic 
spiral, the northern cool part of the dark cloud complex extending between the galaxies, its warm southern 
region, its southernmost star-forming region deficient in radio emission, and the highly polarized northeastern 
ridge associated with the base of an unfolding tidal tail. The whole region of the dark cloud complex shows 
a coherent magnetic field structure, probably tracing the line of collision between the arms of merging spirals 
while the total radio emission reveals hidden star formation nests. The southern region is a particularly intense 
merger-triggered starburst. Highly tangled magnetic fields reach there strengths of ~ 30 /iG, even larger than in 
both individual galaxies, possibly due to compression of the original fields pulled out from the parent disks. In 
the northeastern ridge, away from star-forming regions, the magnetic field is highly coherent with a strong regular 
component of 10 /xG tracing gas shearing motions along the tidal tail. We find no signs of field compression 
by infalling gas there. The radio spectrum is much steeper in this region indicating aging of the CR electron 
population as they move away from their sources in star-forming regions. Modelling Faraday rotation data show 
that we deal with a three-dimensionally curved structure of magnetic fields, becoming almost parallel to the sky 
plane in the southeastern part of the ridge. 

Key words. Galaxies: general - Galaxies; magnetic fields - Galaxies: interactions - Galaxies; individual; Antennae 
- Radio continuum; galaxies 



1. Introduction 

Violently disrupted galaxies show strong departures of gas 
flow patterns from axial symmetry and from symmetric 
spiral shape. Encounters of gaseous disks at velocities in 
excess of 100 km/s give rise to galaxy-scale shocks and 
velocity gradients. In this case magnetic fields can pro- 
vide useful information on the gas flows in the sky plane, 
complementary to any emission-line radial velocity stud- 
ies. In particular the magnetic field can trace regions of 
gas compression (like in Virgo Cluster spirals, Chyzy et al. 
120001 12002|l and of shearing motions stretching the mag- 
netic fields. It is not yet known which gas phases have the 
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strongest influence on the magnetic field evolution. While 
in normal spirals various gas phases seem to be coupled 
via density waves (cold gas in dust lanes, warm and hot 
ionized gas in star-forming regions in spiral arms) in vio- 
lently interacting galaxies various gas phases appear often 
separated. Cold, warm and hot gas may occupy different 
disk regions. This offers an opportunity of studying the 
association of various magnetic field structures with par- 
ticular phases of the interstellar medium. 

The degree to which the field structure follows the gas 
fiows is unknown. While concepts of magnetic field gener- 
ation (see e.g. Beck et al. .1996. Widrow .20021) postulate 
the turbulent magnetic diffusion coefficient to be by three 
orders of magnitude greater than that resulting from flux 
freezing, there are also examples of magnetic fields follow- 
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ing to a large extent the large-scale gas flows (Beck et al. 
I1999|l . In a strongly turbulent environment the magnetic 
field may even become insensitive to strong density wave 
effects (Soida et al. 120011 ■ 

Violently interacting galaxies also constitute a good 
tool to study the global magnetic field evolution. During 
their close encounters the timescale of changing their mag- 
netic field structures is of order of 10* years, much shorter 
than e.g. the timescale of classical dynamo process (e.g. 
Brandenburg & Urpin I1998|l but comparable to that of 
fast dynamos (Moss et al. ll999l Hanasz et al. l2U02|l . There 
is a question to which degree external influences can am- 
plify magnetic fields via the effects of compression and/or 
by providing a strong non-azimuthal magnetic fields in a 
manner similar to the dynamo process. To distinguish be- 
tween the compressed, highly anisotropic fields and unidi- 
rectional ones, generated by the dynamo process, Faraday 
rotation information, thus polarization studies at least at 
two frequencies are needed. 

In this work we present a high-resolution, three- 
frequency radio polarization study of the pair of merg- 
ing spirals NGC 4038/NGC 4039, called "the Antennae". 
This is the first study dedicated to magnetic fields in a 
system of merging galaxies. The galaxies are known to be 
strongly tidally interacting (Toomre & Toomre [l972|l and 
forming two impressive tidal tails (see cartoon sketch in 
Fig-CJ. The neutral gas makes a long extension along the 
southern tail (Hibbard et al. l2000|l but the northern exten- 
sion is devoid of neutral gas. The galaxies are subject to 
a burst of star formation (Liang et al. l2()01l and references 
therein). The young star clusters are distributed along a 
peculiar structure resembling an "inverted 9" (Fig.^, en- 
circling the northern galaxy from the west, then running 
along the eastern edge of the whole system, then turn- 
ing west to joining the southern member of the pair. A 
huge complex of cold dusty matter is visible in the sub- 
mm observations with SCUBA (Haas et al. I2000|l in the 
"overlap area" of the galaxy disks concentrated around 
RA2000 = 12'^01'"54!7 Dec2ooo = -18°52'52" (denoted by 
a cross in Fig.^. It coincides also with a strong CO peak 
found by Stanford et al. The mid-infrared ISO 

observations (Mirabel et al.|1998') shows that the most in- 
tense starburst occurs in optically obscured place to the 
east from the southern galaxy where emission comes from 
gas and dust heated by massive stars. The ionized gas 
in NGC 4038/39 shows a complex velocity field (Amram 
et al. 11992(1 while the X-ray emitting hot gas has a pe- 
culiar distribution with the evidence of coUimated out- 
flows (Read et al. .1995 ) . High resolution Chandra ACIS 
data (Fabbiano et al. 12001(1 reveals a population of ex- 
tremely luminous point-like sources (probably black hole 
binaries) and soft thermal X-ray emission associated with 
star-forming knots but sometimes also intermingle with 
warm (Ha) gas. 

The galaxy pair NGC 4038/39 has been studied using 
VLA at 1.47 GHz and 4.89 GHz by Hummel & van der 
Hulst fl986') . They found a radio ridge extending along the 
mentioned "inverted 9" structure, however, their maps do 
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Fig. 1. Sketch of the NGC 4038/39 merging system with 
an enhanced DSS optical image in the background. The 
black cross marks a concentration of large complex of dark 
dusty clouds in the "overlapping region" . 



not show much extended structure, moreover the authors 
did not analyze the polarization information. 

In this paper we present maps of NGC 4038/39 in the 
VLA D-configuration at 8.44 GHz, combined C-|-D-array 
maps at 4.86 GHz as well as C-array data at 1.49 GHz, 
much more sensitive to extended diffuse emission than the 
earlier maps by Hummel & van der Hulst (1986). After 
presenting some observational aspects of our radio data 
in Sect. [2 we introduce the total radio and polarized 
intensity maps in Sect. |2I together with the distribution 
of spectral index, Faraday rotation and depolarization. 
Decomposition of thermal and nonthermal radio intensity 
as well as polarization data results in identifying distinct 
radio-emitting regions (Sect.QJ. By means of profiles along 
different direction across the merging system we study the 
association of magnetic fields with different gas phases us- 
ing available data in radio continuum, H i, optical. Ha, 
X-ray, and dust emission fSect. UTTjl . In Sect. 14.21 we dis- 
cuss total and regular magnetic field strengths derived in 
various parts of the interacting system and compare them 
to values known from other galaxies. We give some idea 
of the magnetic field origin and its regularity in Sect. 14.31 
Finally, we present a tentative scenario of the gas and 
magnetic field evolution in the Antennae in Sect. 14.41 The 
summary and suggestions of further works are given in 
Sect. El 

2. Observations and data reduction 

Observations were made at 8.44 GHz and at 4.86 GHz, 
using the VLA of the National Radio Astronomy 
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Fig. 2. Total intensity contours of radio emission of 
NGC 4038/39 at 4.86 GHz made from a combination of 
VLA C and D-array data overlaid upon the optical im- 
age from HST (courtesy B.C. Whitmore from STScI). 
The contour levels are 0.2, 0.4, 0.7, 1.1, 1.7, 2.9, 4.5, 7.4 
mjy/b.a. The rms noise is 0.021 niJy/b.a. and the beam 
size is 6'.' 5 x 4'.' 5 (uniform weighting). White circles mark 
the positions of galaxies cores and the black cross a con- 
centration of cold dust in the "overlapping region" . 

Observatory-'^ in the compact D-array configuration. 
Additional observations at 4.86 GHz and the measure- 
ments at 1.49 GHz were performed with the C-array. The 
observation times in C-array were 5.5 hours at 1.49 GHz 
and 6 hours at 4.86 GHz, while in D-array the observation 
times were 5 hours at 8.44 GHz and 5.5 hours at 4.86 GHz. 

The intensity scale at all frequencies was calibrated by 
observing 3C286, adopting the flux densities of Baars et al. 
P977|l . The position angle of the linearly polarized inten- 
sity was calibrated using the same source with an assumed 
position angle of 33°. At 4.86 GHz and 8.44 GHz the 
calibrator 1157-215 was used to determine the telescope 
phases and the instrumental polarization. To check these 
calibration procedures 3C138 was observed once each day. 
At 1.49 GHz the phase calibrator was 1156-221. 

The data reduction has been performed using the AIPS 
data reduction package. The edited visibility data, cali- 
brated and self-calibrated in phase, were Fourier trans- 
formed to obtain maps in Stokes parameters I, Q and U 
at three frequencies. At 4.86 GHz the maps were made us- 
ing a combination of C and D-arrays with various weight- 
ings for particular purposes. The combined data were self- 
calibrated in phase. To increase the sensitivity to extended 
structures we merged the I, Q and U data at 8.44 GHz 

^ NRAO is a facility of National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



in the UV plane with observations at 10.45 GHz made 
by us with the 100-m Effelsberg radio telescope^. The 
brightness values at 10.45 GHz were rescaled to 8.44 GHz 
assuming a spectral index of 0.7. Faraday rotation mea- 
sures only scarcely exceed lOOradm"^ (Fig-EJ which cor- 
responds to a polarization angle offset of 3° between 10.45 
and 8.44 GHz. By comparing the merged VLA/Effelsberg 
maps and pure VLA ones we state that any possible po- 
larization angle offset between 10.45 GHz and 8.44 GHz 
does not introduce visible artifacts in the merging process. 

Finally the Q and U maps at three frequencies were 
combined to get distributions of the linearly polarized in- 
tensity (corrected for the positive zero level offset) and 
of the position angle of the apparent magnetic vectors 
(B-vectors; not corrected for Faraday rotation observed 
E- vectors rotated by 90°). 



3. Results 

3.1. Total radio intensity 

Details of the distribution of total radio intensity in 
NGC 4038/39 are best visible in our combined high- 
resolution C-f D-array map at 4.86 GHz (Fig. EJ. The 
bright emission in the NW region follows the "inverted 9" 
feature as already described by Hummel & van der Hulst 
(1986). Local total intensity peaks coincide with the nu- 
clear regions of both galaxies. An extended, radio-bright 
feature is located in the southern part of mentioned large 
complex of dark clouds where a mixture of dark patches 
and Hii regions is present (Whitmore et al. I1999|l . While 
localized peaks along the "inverted 9" in the NW por- 
tion coincide with particular star-forming complexes, the 
mentioned radio-bright region does not show particularly 
bright Ha-emission and its peaks seem to be loosely as- 
sociated with optically bright regions (Fig. |2J). No radio 
emission was detected from the gas-poor northern tidal 
plume. 

The low-resolution total intensity maps of 
NGC 4038/39 at 8.44 GHz and 4.86 GHz (Fig. ^ 
show consistently a strong emission peak in the southern 
part of the dark cloud complex extending between the 
galaxies. Some weak maxima coincident with their nuclei 
are seen as well. Both maps show also an extension 
along the southern gas-rich tidal tail. No emission is 
detected from the gas-poor northern tail. A weak tongue 
of emission towards the southwest is seen consistently 
at both frequencies and is probably real. The total 
intensity maps at both frequencies also possess a steep 
brightness gradient along the NE straight boundary of 
the interacting pair. 



^ The 100-m telescope at Effelsberg is operated by the Max- 
Planck-Institut fiir Radioastronomie (MPIfR) on behalf of the 
Majc-Planck-Gesellschaft 
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Fig. 3. Left panel: the total radio intensity contours and apparent B-vcctors (not corrected for Faraday rotation) 
of polarized intensity of NGC 4038/39 at 8.44 GHz made from a combination of VLA and Effelsberg data with a 
resolution of 17" x 14" (natural weighting), overlaid upon a DSS image. The contour levels are 0.07, 0.11, 0.33, 0.77, 
1.8, 3.3, 6.6, 11, 18 mJy/b.a., the rms noise is 0.022 mJy/b.a., a vector of length of 100" corresponds to a polarized 
intensity of 1.25 mJy/b.a. Right panel: the contours of total intensity and apparent B-vectors of polarized intensity 
of NGC 4038/39 at 4.86 GHz made from a combination of VLA data in C and D with a resolution of 17" x 14" 
(naturally weighed and convolved) overlaid upon the DSS image. The contour levels are 0.05, 0.12, 0.30, 0.53, 1.2, 2.1, 
3.3, 5.3, 9.0, 17, 24 mJy/b.a. and the rms noise is 0.015 mJy/b.a., a vector of length of 100" corresponds to a polarized 
intensity of 1.25 mJy/b.a. 



3.2. Polarization maps 

The maps of polarized intensity at 8.44 GHz and 4.86 GHz 
(Fig. are suggestive for several structural components. 
In the NW disk ( "inverted 9" ) the polarized emission co- 
incides with the optical one (Fig.Oright panel). However, 
the orientations of apparent B-vectors do not follow the 
"inverted 9" feature. Instead, they form a rudimentary spi- 
ral shifted westwards from the northern galaxy nucleus, 
being centered on RA= 12'^1™51!6 DEC= -18°51'57" 
(Fig. |3| right panel, a white cross). The B-vectors at both 
frequencies are inclined by some 45° with respect to the 
total intensity ridge, thus, this feature is not caused by 
Faraday rotation effects. Close to the centre of north- 
ern galaxy the B-vectors are almost radial, following local 
dust lanes visible on HST image in Fig. [21 Some traces of 
another spiral with a very large pitch angle are present 
around the nucleus of the southern galaxy. They are best 
visible at 8.44 GHz (Fig. El left panel). 

In the NE region of the interacting system the apparent 
B-vectors run parallel to a rather sharp, linear NE bound- 
ary of the optical emission. The map of polarized intensity 
at 8.44 GHz, much less affected by Faraday effects, shows a 
bright peak at RA= 12i^l™55!5 DEC= -18°51'57" (Fig.El 
left panel, black cross). The polarization degree increases 



along this feature from about 5% in its northern part to 
about 40% at its southernmost tip. The polarized inten- 
sity at 4.86 GHz shows more clearly the southern part of 
this structure and forms a bright ridge extending along the 
NE boundary of the optical image of the system (Fig. 01 
right panel). 

Another polarization component apparently extends 
along the complex of dark clouds visible in the B /K colour 
index image (Hibbard et al. i2000^) shown in Fig. 01 (left 
panel). It is separated from the polarized ridge discussed 
above by a shallow depression, weakly visible at 8.44 GHz, 
but becoming much more conspicuous at 4.86 GHz. The 
northern part of the dark dust complex, devoid of observ- 
able Ha emission (Whitmore et al. I1999|l coincides with a 
local maximum of polarized intensity (see Fig.^left panel, 
white cross). However, at this position the polarization 
degree is rather moderate: about 10% at 8.44 GHz and 
only 6% at 4.85 GHz. The polarized emission at 8.44 GHz 
extends southwards along the region of strong reddening 
towards the bright total intensity peak (marked in the fig- 
ure by a white circle). The polarization degree decreases 
southwards, becoming only 1% at 8.44 GHz at the position 
of the total intensity maximum. The polarized intensity at 
4.86 GHz decreases southwards even faster. 



K. Chyzy et al.; Magnetic fields in the Antennae 



5 




12 0202 00 01S8 56 54 52 50 48 1 2 02 02 01 58 56 54 52 50 48 

RIGHT ASCEHSlOn (J2000) RIGHT ASCCHSIOH (J2000) 



Fig. 4. Left panel: The contours of polarized intensity with apparent B-vectors of polarization degree of NGC 4038/39 
at 8.44 GHz overlaid upon the B/K colour index map (where dark areas indicates blue colour) from Hibbard et al. 
1)2000(1 . The contour levels are 0.04, 0.13, 0.26, 0.46 mJy/b.a. and the rms noise is 0.013 mJy/b.a. polarization degree 
of 12.5%. The beam size is 17" x 14". Right panel: contours of polarized intensity of NGC 4038/39 and apparent 
B-vectors of polarization degree at 4.86 GHz overlaid upon the DSS image. The contour levels are 0.04, 0.11, 0.21, 
0.28, 0.42 mJy/b.a and the rms noise is 0.014 mJy/b.a. A vector of 10" corresponds to a polarization degree of 12.5%. 
The beam size is 17" x 14". 
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Fig. 5. The distribution of spectral index obtained by 
fitting a power low to the values of total intensity at 8.44, 
4.86 and 1.49 GHz in NGC 4038/39. The maps of total 
intensity were convolved to a common beam of 22" x 14". 
The contours are 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1. 



3.3. Spectral index distribution 

The spectral index distribution in NGC 4038/39 (Fig. 01 
has been obtained by fitting for each map pixel a power 
law oc to the values of total intensity at three fre- 
quencies. The discussed interacting pair is characterized 
by strong variations of the spectral index over its surface. 
The western part of the "inverted 9" has a rather flat spec- 
trum with a ~ 0.6, flattening northwards to 0.53. The ra- 
dio spectrum also gradually flattens along the dark cloud 
complex in the central region, changing from about 0.67 
at its northern end to 0.58 in the radio bright region co- 
incident with the southern part of dust complex. Nuclear 
regions of both galaxies also have a rather flat spectrum 
with a of about 0.64. 

In contrast to that the gap between the galaxies, show- 
ing also depression in total radio intensity (Fig. |3Jl has a 
spectral index of about 0.94. A very steep spectrum has 
been found along the northeastern polarized region: most 
of it has a > 1, steepening towards the tidal tail to values 
of 1.2-1.3. There is also some spectral steepening across 
the ridge: from a ~ 0.90-0.93 at its interface with a central 
cloud complex to values above 1.0 at the outer boundary. 

3.4. Distribution of Faraday rotation and 
depolarization 

Maps of Faraday rotation and depolarization in 
NGC 4038/39 (Figs. El and 0) have been computed be- 
tween the frequencies of 4.86 and 8.44 GHz using the 
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naturally weighted maps at the resolution of 14" x 17". 
The galaxies lie at 287° Galactic longitude and 42.4° 
Galactic latitude where Simard-Normandin & Kronberg 
H198U|I show low foreground rotation measures changing 
sign from place to place and absolute values < 30 radm~^. 
For this reason no correction for the foreground Faraday 
rotation has been applied. 

Over most of the eastern part of NGC 4038/39 and 
in the westernmost part of "inverted 9" the Faraday ro- 
tation measures (RM) are generally small. Their absolute 
values in these regions exceed 50 rad m^ only in small ar- 
eas comparable to the beam size (Fig. 01 . The dividing 
line between the positive and negative values may depend 
on possible bias due to the foreground Faraday rotation. 
However, some statements do not depend on the uncer- 
tainty of foreground RM by some ±30radm~^. We note 
first the vertical belt of large negative RM's (between —50 
and — 150radm~^), about 30" wide, running north-south 
through the nucleus of the northern galaxy. Further to the 
south, around the southern galaxy nucleus, sudden RM 
jumps to strongly positive values up to > -1-100 rad m~^ 
are visible. We also note a substantial coherence of the sign 
of RM values. Independently of the assumed foreground 
rotation we get invariably a large domains of several beam 
areas having a constant sign of Faraday rotation measure 
and hence of the magnetic field direction. 

In the eastern region of NGC 4038/39 the Faraday 
depolarization (FD=polarized intensity at 4.86GIIz / po- 
larized intensity at 8.44GIIz; Fig. forms a valley with 
FD ~0.2-0.3 along the channel occupied by the hot gas 
detected there by Chandra (Fabbiano et al. I2()()l(l . Very 
strongly depolarized regions (FD< 0.2) also appear near 
the nuclei of both galaxies. Another strong Faraday de- 
polarization (FD~ 0.12) occurs around RA= 12'>01'°54!0 
DEC= -18°51'35" (marked by a black cross in Fig. □) 
where the Faraday rotation measure changes locally from 
-120radm-2 to -h50radm-^ 



4. Discussion 

To study the physical conditions and evolutionary stages 
of various locations in the Antennae system and interrela- 
tions of magnetic fields with particular phases of the inter- 
stellar medium, a comparison of our radio data with other 
spectral domains is necessary. For this purpose, mainly 
based on our radio maps in total and polarized intensity, 
we define several distinct regions labelled by letters A-F in 
Fig. |S1 Region A is located in the middle of the polarized 
northeastern ridge, at the point where a bright ridge of po- 
larized intensity coexists with a similarly oriented feature 
in Hi fFig. 113(1 . Regions B, C and D denote various parts 
of the so-called "overlap region" , constituting a dark cloud 
complex extending between the two interacting galaxies. 
Region B lies at its northern end while regions C and D 
denote the radio total intensity peaks at its southern tip 
(Fig. 121). As region E we define the middle of the extended 
pool of X-ray emitting gas (Fig. El Fig- IT^ . As region F 
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Fig. 6. Faraday rotation measure distribution in 
NGC 4038/39 computed between 8.44 GHz and 4.86 GHz 
plotted as symbols (see the Figure legend) overlaid upon 
the grayscale plot of rotation measure errors. The map 
resolution is 17" x 14". Black crosses mark the positions 
of the galaxy cores. 



we term the NW part of the relic spiral arm in the NW 
region of the system. 

In the following subsection, using our total intensity 
maps at 8.44, 4.86 and 1.49 GHz, we first construct 
and present maps of thermal and nonthermal intensity 
of NGC 4038/39. We then analyze in next subsections the 
local interrelations of different ISM phases in regions A-F 
by comparing various distributions of radiation of different 
species along and across these particular regions. 

4.1. Thermal and nonthermal intensity in 
NGC 4038/39 

Our total intensity maps at 8.44, 4.86 and 1.49 GHz allow 
to separate the thermal and nonthermal radio emission 
components, assuming a constant nonthermal spectral in- 
dex ant over the whole radiation area. However, there are 
local exceptions from this assumption where spectral in- 
dex become particularly steep. There are regions in which 
the observed spectral index reaches the value of 1.0 (the 
radio "valley" between the galaxies on the western side 
of the system) or even 1.2-1.3 (the base of the tidal tail. 
Fig. EJ, probably because of strong energy losses by CR 
electrons. 
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Fig. 7. Faraday depolarization distribution in 
NGC 4038/39 computed between 8.44 GHz and 4.86 GHz. 
The contours are 0.2, 0.4, 0.6, 0.8. The map resolution is 
17" X 14". 
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Fig. 8. Total radio intensity map and apparent B-vectors 
of polarized intensity of NGC 4038/39 at 4.86 GHz made 
from a combination of C and D array data, showing the 
location of particular radio emission regions (A-F) as dis- 
cussed in Sect.0] The map resolution is 12" x 6". Shadow 
areas denote slices along which the profiles of various 
species were obtained (see Sect.^: a) a slice along north- 
eastern ridge; b) a slice across the ridge; c) a slice along 
the "overlapping region" . 



To compute the maps of thermal and nonthermal in- 
tensity (and thus also the map of thermal fraction) we 
tried to find the best value which could be adopted for 
mean ant in possibly largest area of the disks ignoring 
mentioned regions of strong electron energy losses. For this 
purpose we used the integrated spectrum of NGC 4038/39 
(Fig.|5J) constructed from the available literature data and 
our own measurements. The flux density at 10.45 GHz 
was obtained by integrating our Effelsberg map in circu- 
lar rings. The spectrum can be equally well represented 
by a whole range of combinations of nonthermal spectral 
index and mean thermal fraction at 10.45 GHz ranging 
from ant = 0.69 fth = to a^t — 0.8 and fth — 0.5. 
All fits between these values lie well within 3a flux den- 
sity errors (Fig. Ej) and yield very similar r.m.s. residuals. 
A further increase of a mean thermal fraction (hence a 
steeper nonthermal spectrum) causes a rapid increase of 
residuals and yields an unrealistically curved total spec- 
trum. On the other hand, adopting a flatter nonthermal 
spectrum (ant < 0.8) seems unreasonable in computing 
the maps because large disks areas are occupied by a spec- 
trum steeper than this value. We used thus a compromise 
assumption of ant=0.8 and computed maps of thermal 
and nonthermal intensity as well as the map of thermal 
fraction at 8.44 GHz, excluding from the analysis regions 
with a steeper spectrum. While this assumption is valid 
for the area covered by optically bright galaxies, the non- 



thermal spectrum is definitely steeper in the vicinity of 
the radio extension along the northeastern tail and in the 
radio "valley" between the galaxies. 

We note here that a mean thermal fraction of 0.5 at 
10.45 GHz implied by this assumption is larger than the 
mean value of 0.30±0.05 at 10 GHz estimated for 74 galax- 
ies from the Shapley- Ames survey by Niklas et al. H1996|l . 
This is not surprising as the pair NGC 4038/39 actively 
forms stars and produces ionized gas very efhciently. 

Regions of substantial thermal fraction (ff^ at 
8.44 GHz of 30%-40%) are found at the positions of the 
nuclei, in the western part of the "inverted 9" (region 
F), and along the dark cloud complex. In the southern 
part of the obscured regions C and D the thermal fraction 
reaches some 47%. The distribution of thermal brightness 
at 8.44 GHz is compared to the distribution of H ii regions 
in Fig.^l The ridge of enhanced thermal intensity follows 
well the chain of H ii regions in the "inverted 9" , making 
some minimum east of the northern galaxy nucleus, (close 
to our region B) then following the belt of Hii regions 
along the central dark cloud. The thermal brightness has 
a maximum in the southern part of the cloud, while no 
particular peak in Ha emission nor any concentration of 
Hii regions is seen. This is discussed further in Sect. 14.131 

The nonthermal intensity distribution (Fig. 1111) is 
smoother and fills the whole system, as expected from 
smearing small-scale features by the cosmic ray propaga- 
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NGC 4038^9 Nonthermal emission + H alpha 



Fig. 9. The integrated spectrum of NGC 4038/39 com- 
piled from the hterature data: 1 - Wright & Otrupcek 
(CnHjl, 2 - Wright & Otrupcek (TM]|, 3 - Hummel & 
van der H ulst ijWgBa . 4 - Condon ( 19^ . 5 - Wright & 
Otrupcek fT^. 6 - Griffith et al. lITTMIl . 7 - Wright & 
Otrupcek H1990|l . 8 - this paper. The dashed line shows the 
model spectrum assuming a thermal fraction at 10.45 GHz 
of 50% and a nonthermal spectral index of 0.8. 
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Fig. 10. The contour map of the thermal intensity at 
8.44 GHz in NGC 4038/39 overlaid upon the grayscale im- 
age in the Ha line made by HST (courtesy B.C. Whitmore 
from STScI). The contour levels are 0.1, 0.2, 0.7, 1.3, 2.2, 
6.6, 8.8 mJy/b.a. and the map resolution is 22" x 14". 



tion. The southern part of the central dark cloud shows 
also a substantial enhancement of nonthermal intensity. 

4.1.1. Star formation in the vicinity of the 
northeastern ridge 

The profile of the total and polarized radio intensity com- 
pared to that of Hi, X-ray and the radio thermal and 
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Fig. 11. The contour map of nonthermal intensity at 
8.44 GHz from NGC 4038/39 overlaid upon the grayscale 
image in the Ha line made by HST (courtesy B.C. 
Whitmore from STScI). The contour levels are 0.1, 0.2, 
0.7, 1.3, 2.2, 6.6, 8.8 mJy/b.a. and the map resolution is 
22" X 14". 



nonthermal intensity across the northeastern ridge is pre- 
sented in Fig. El Most of species peak at the distance from 
the slice origin x ~ 80" (region B in Fig. ^ . This corre- 
sponds to the northern part of the dark cloud complex 
extending between the galaxies, flanked on its NE side by 
a star-forming region. The separation of thermal and non- 
thermal intensity is visible in the middle panel of Fig. 1141 
The high resolution optical image (shown in Fig. [21) helped 
us to associate the thermal emission peak with the star- 
forming clump and the nonthermal one with a dusty cloud. 
This means that there is not much hidden star formation 
in the cloud itself. Instead there is a strong, partly reg- 
ular magnetic field (note the peak of polarized intensity) 
into which the electrons produced in the star-forming re- 
gion propagate, yielding enhanced nonthermal intensity. 
We note that the distribution of the A15^m brightness 
(obtained by Mirabel et al. 119981 with ISO) agrees with 
the total (i.e. sum of thermal and nonthermal) intensity, 
its peak being located between the maxima of these two 
emission components. Both thermal and nonthermal in- 
tensity decrease towards the NE. 

The X-ray emission peaks about x — 100" in a pool of 
hot gas (region E) . This "X-ray ridge" does not correspond 
to any increase in total radio intensity. According to Zhang 
et al. I|2001|l this region is filled mostly with intermediate 
age (3 < i < 16 Myr) stellar clusters. It may be a relic of 
an old star formation burst where the cosmic ray electrons 
lost most of their energy, making no significant contribu- 
tion to the nonthermal intensity. Synchrotron lifetime of 
electrons with ~ 4 GeV energy in a 23 /l(G field (Table ^1 
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Fig. 12. Contour map of polarized intensity and apparent 
B-vectors of polarization degree at 4.86 GHz overlaid upon 
the grayscale image of NGC 4038/39 in X-rays observed 
by Chandra satellite (Fabbiano et al. I2001|l . The contour 
levels are 0.04, 0.11, 0.21, 0.28, 0.42 mJy/b.a. and the map 
resolution is 17" x 14". 

is ~ 5 X 10^ yr which gives a lower limit for the age of the 
star burst. 

Polarized intensity in region E shows a striking min- 
imum at 8.44 GHz and at 4.86 GHz (see Fig. ^ to- 
gether with strong wavelength-dependent depolarization 
(FD "valley" in Fig.d). At 4.86 GHz the degree of polar- 
ization is reduced by a factor of 0.35 relative to 8.44 GHz. 
As Faraday rotation there is smaller than 50 rad m~^ 
(Fig. 5), depolarization due to differential Faraday rota- 
tion cannot be more than by a factor of 0.98. Thus, the 
wavelength-dependent depolarization is due to Faraday 
dispersion CBurn [TUm Sokoloff et al. ITOMjl . A RM dis- 
persion of CT^j^ of 1.1 X 10^ yields a depolarization of 0.30 
at 4.86 GHz and 0.84 at 8.44 GHz giving the observed ratio 
of 0.35, but insufficient to explain the minimum in polar- 
ized intensity in region E at 8.44 GHz. In that case ad- 
ditional wavelength-independent beam depolarization by 
an increased degree of small-scale magnetic field tangling 
is required. 
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Fig. 13. Upper panel: contours of polarized intensity 
at 4.86 GHz in VLA C-array with uniform weighting 
(showing the brightest polarized structures at highest res- 
olution) overlaid upon to the distribution of Ha emis- 
sion (grayscale). The cross-hatched area shows the region 
where HI emission was observed by Hibbard et al. I)2000|l 
at the resolution of 11" x 7"4 exceeds 100 mJy/b.a. Lower 
panel: the contours of polarized intensity and grayscale of 
Ha emission (as above) and the region occupied by the dif- 
fuse X-ray emission (hatched area, Fabbiano et al. l2()()l|l . 



4.1.2. Magnetic fields and cosmic rays in the 
northeastern ridge 

At a; ~ 120" (region A) a peak of polarized intensity is 
observed, coinciding with a H i maximum (Fig. 1141 see 
also Fig. 113(1 . It does not correspond to any rise of other 
constituents. According to recent numerical simulations 
kindly supplied to us by Dr Englmaier (priv. comm., see 
also Englmaier et al. I2001|l this region constitutes the be- 
ginning of the tidal (both stellar and gaseous) tail, un- 
folding towards the east and northeast. Because of little 



star formation this region does not correspond to any par- 
ticular total intensity features, instead it is well revealed 
by the polarized intensity because of highly ordered mag- 
netic fields. According to the suggestion by Fabbiano et 
al. ipnHjl this region may be subject to external compres- 
sion by the ambient gas pushing inwards the pool of hot 
gas. This may imply that our region A, with well aligned 
polarization vectors in front of the X-ray maximum, could 
be the result of magnetic field compression. Gas infall 
in this region has also been suggested by Hibbard et al. 
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NGC 4038/39 profiles across the ridge 



NGC 4038/39 profiles along the ridge 
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Fig. 14. The profiles of various species across the NE 
nonthermal ridge integrated along the slice b) marked in 
Fig.lHl 

H2UUU|I . However, the high resolution map of polarized 
intensity (Fig. I13|) does not show any narrow ridge, as 
expected in case of external compression. Moreover, the 
highly polarized region is significantly displaced southeast- 
wards with respect to the X-ray ridge and also to the H i 
ridge fFig. I13|l . It is more likely that the increased degree 
of magnetic field regularity is caused by field geometry and 
decaying turbulence and/or shearing gas motions (see be- 
low) than by gas infall pushing inwards the X-ray emitting 
matter. 

Fig. El shows the distribution of various constituents 
along the discussed polarized ridge. There is still some 
weak star formation in its N part, as indicated by a total 
intensity peak at a; ~ 80". The polarized intensity peaks 
more southwards at ~ 100" at 8.44 GHz (Fig. [Tsl) and 
a; ~ 120" at 4.85 GHz. While both total and polarized in- 
tensity decrease towards the tidal tail, the polarized inten- 
sity fades somewhat farther away than the total intensity 
(Fig. El lower panel). This means that the magnetic field 
gets more ordered (the polarization degree is increasing) 
as the turbulent motions decay away from star-forming re- 
gions. A substantial part of a high polarization of this re- 
gion may come from a "pseudo-regular field" , i.e. a highly 
anisotropic random one, stretched by shear accompanying 
the unfolding motion of the tidal tail. Such a field struc- 
ture yields a high degree of polarization, not accompanied 
by strong Faraday rotation, because the stretching of ran- 
dom field alone cannot generate coherent uni-directional 



Fig. 15. The profiles of various species along the NE 
nonthermal ridge integrated along the slice a) marked in 
Fig.lHl 



magnetic fields. We found indeed that Faraday rotation 
in the southern, weakly star-forming part of the polarized 
ridge is small (< 30radm~^, see Fig-jSJ. 

However, the Faraday rotation measures (in absolute 
values) increase above lOOradm^^ along the ridge of po- 
larized emission towards the north. These strongly non- 
zero Faraday rotation indicate that we deal with gen- 
uinely regular magnetic fields. Together with the increase 
in RM, the degree of polarization decreases. This anticor- 
relation suggests that we deal with the same large-scale 
coherent magnetic field and variations in the orientations 
of magnetic lines are responsible for the observed vari- 
ations of RM and the field regularity. This explains the 
high Faraday rotation in the north. The small degree of 
polarization (low field regularity) in the north is in line 
with the high star formation rate, a relation also observed 
in other regions of NGC 4038/39 (see Sect. ETSj). We 
believe that the southern region also contains strong co- 
herent fields but they run almost parallel to the sky giving 
little Faraday rotation. 

A simple model reproducing simultaneously Faraday 
rotation measures, polarized intensities and thermal in- 
tensity at 8.44 GHz has been constructed. We assume 
that 20% of the observed thermal intensity (Fig. I10|l is 
due to diffuse ionized gas responsible for Faraday rotation. 
Varying the pathlength of thermal emission between 100 
and 400 pc, the nonthermal pathlength between 1 kpc and 
2kpc, and the ionized gas filling factor between 0.05 and 
0.1 (see Ehle & Beck ll955)l we found that all radio data 
in the SE part of the polarized ridge (marked by a black 
cross in Fig. (HJ require a regular field of about 10 ± 1 /xG 
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(see TableP), inclined to the sky plane by some 4±2°, only 
weakly depending on the assumptions listed above. This 
result suggests that most of the regular magnetic field lies 
in the sky plane. 

In the northern part of the ridge (marked by a black 
circle in Fig. |Hl, showing higher RM, the same analysis 
implies that the regular field of about 9 ± 2 fiG is running 
away from the observer at the angle of 59 ± 9° to the 
sky plane. From this place towards the SE (in a region 
marked by a black triangle in Fig. ISJ the regular magnetic 
field changes its direction. It remains inclined to the sky by 
some 34 ± 9° but points towards us. Our model indicates 
that we generally deal with a three-dimensionally curved 
structure of magnetic fields. 

These results suggests that we need a higher degree of 
coherence in the southern part of the NE ridge as expected 
for decaying turbulent motions away from star-forming 
regions. While the above mechanism would be the main 
cause of coherent fields, the shearing motions may align 
the regular field along the tail. 

Along the nonthermal ridge the radio spectrum 
(Fig. I15|l gradually steepens (also the nonthermal one as 
discussed in Sect. I4.ljl . indicating an aging population of 
CR electrons as they propagate further away from their 
sources in star-forming regions. The proximity of a hid- 
den strong star formation nest in the southern part of the 
"overlapping region" (D in Fig. IHJ does not result in an 
increased CR electron content in the ridge, probably be- 
cause of a strong regular magnetic field perpendicular to 
the required direction of CR propagation. 

4.1.3. The dark cloud complex 

Fig. El shows the brightness variations of various con- 
stituents along the dark cloud complex extending to SE 
from the nucleus of the northern galaxy (seen on e.g. 
the B/K colour index map in Fig.01 left panel). Starting 
from the northwest there is a small peak in polarization 
(x ~50"), corresponding to the interface between the rudi- 
mentary magnetic spiral in the NW disk (see Sect. l3.2(l and 
the polarized ridge discussed in Sect. 14.1.11 The whole 
relic spiral associated with the "inverted 9" optical fea- 
ture deserves attention, too. In this region the observed 
B-vectors have almost a constant inclination (~ 30°-40°) 
with respect to the narrow blue optical ridge, thus the 
magnetic field cannot be simply aligned by compression 
as this structure unfolds. We suspect this feature to be 
the relic of the galaxy's original magnetic field. However, 
Faraday rotation measure in the area of this rudimentary 
spiral undergoes strong variations and jumps, changing 
suddenly between ~ — lOOradm^^ and ~ -|-100radm^^. 
This indicates that the magnetic spiral has in fact a com- 
plex three-dimensional structure as does also the velocity 
field fBarnes fl988|l . Similar but weakly visible spiral with 
a large pitch angle is visible also around the southern nu- 
cleus. 
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NGC 4038/39 profiles along the overlaping region 
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Fig. 16. The profiles of various species in the "overlap- 
ping region" integrated along the shce c) marked in Fig.|Hl 



Moving to the southeast, at x ~ 75" fFig. Ilt)|) there is 
a small peak in almost all species (except polarization and 
H i) corresponding to an enhanced star formation around 
the nucleus of the northern galaxy. The polarized intensity 
continues to rise towards the SE, reaching a strong peak 
(as does the Hi) at x c^i 85" (region B) while all other 
quantities, sensitive to the star-forming activity show a 
minimum close to this place. This region corresponds to 
the northern part of the dark cloud complex. Our observa- 
tions imply thus that there is no hidden strong star forma- 
tion there. In fact the sub-mm observations with SCUBA 
(Haas et al. |^00) show that this region is radiating quite 
strongly at A850/im being very weak at A450/im, which 
indicates its low dust temperature. 

More to the south at x ~ 130" the radio thermal and 
nonthermal brightness, as well as the intensity at A15//m 
(ISO) show a strong peak fFig. I16|l . constituting a beam- 
smeared blend of regions C and D as defined in Fig. |H1 
The peaks of the radio and infrared emission are more 
than twice stronger compared to those at the position of 
the nucleus of northern galaxy. In contrast, the Ha and 
X-ray peaks are considerably weaker than at the northern 
nucleus. In case of Ha emission this means strong extinc- 
tion in the optical domain, hiding a very strong burst of 
star formation. This agrees with results by Haas et al. 
f2000) who find there strong emission at A450/im indicat- 
ing considerably warmer dust. The peak of thermal radio 
intensity is particularly appealing. The radio thermal in- 
tensity can therefore serve as a good tracer of thermal gas 
ionized by recent star formation, but still hidden in the 
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optical domain by heavy extinction in a dusty molecular 
environment. 

Whitmore & Zhang p999|l argued that most radio 
emission from the radio-bright phase of young star clus- 
ters in the Antennae originate from Ha-emitting gas. Our 
results show that their conclusion is not valid for the ex- 
tended radio emission. Indeed, the compact sources com- 
prise only 11% of the total radio intensity studied here. 
The distinct behavior of Ha and total radio intensity is 
visible for example in Fig. 1161 where their profiles along 
the "overlapping region" show different shapes. 

To explain a weak X-ray emission in the southern part 
of the dark cloud complex we have two possibilities: ei- 
ther strong extinction by a H i concentration, or the star 
formation burst may be so young that did it not yet pro- 
duce a large reservoir of diffuse hot gas but sufficiently old 
to produce radio-emitting CR electrons, hence its age is 
between 10^ yr and 10^ yr. The discrimination between 
these possibilities is impossible without analyzing the X- 
ray spectra, which is beyond the scope of this paper. 

The degree of polarization and the polarized intensity 
is highest in the cool northern part of the "overlapping 
region" (close to B). They both decrease gradually to- 
wards the SE while the total (i.e. the sum of thermal and 
nonthermal) intensity continuously rises (Fig. I16|l . This 
implies that in the southern region of the dark cloud com- 
plex the regular magnetic field becomes more efficiently 
destroyed by hidden star formation. We also note that the 
B-vectors keep a coherent structure along the whole re- 
gion (Fig. 01). We deal here with genuinely regular fields: 
the Faraday rotation measure keeps constantly high ab- 
solute values (> lOOradm"^) in two large coherent do- 
mains: negative in the northern part, changing to pos- 
itive in its southern region (Fig. EJ. The analysis done 
in the same way and under the same assumptions as in 
Sect. ^^21 shows that equipartition regular magnetic field 
in the northern part (region B) is about i?reg — 5 fiG (as 
estimated from the polarized brightness at 8.44 GHz) and 
is inclined to the sky plane by some 37±5°. This is implied 
by the Faraday rotation measure and thermal gas density 
derived from the thermal intensity. The sudden change of 
RM in the middle of the cloud complex also means a sign 
reversal of the line-of-sight regular field component. We 
tentatively speculate that this might be a meeting point 
of the global regular magnetic fields of the two galaxies. 

As also noted by Haas et al. I)2U(JU|I the southernmost 
tip of the dark complex (region C, distinguished in our 
high resolution total radio intensity map. Fig. ^ has a 
several times higher ratio of A15/im-to-radio intensity than 
region D. This is also visible in Fig. 1161 fupper panel) as a 
shift of the infrared peak with respect to the radio emis- 
sion. It is possible that our region C contains a hidden, 
strongly star-forming nest residing in the environment of 
weaker total magnetic field or the star formation burst is 
so young (< 10^ yr) that not enough CR electrons have 
been produced so far. 



4.2. Magnetic field strengths 

For the detailed analysis of magnetic fields in 
NGC 4038/39 we performed calculations of total i?tot 
and regular -Bicg equipartition magnetic field strengths 
and of the degree of field regularity i?rcg/Sian (the ratio 
of regular to random field strengths) in different parts 
of the interacting system (Table For all calculations 
we used maps of total and polarized radio intensity 
at 8.44 GHz. We adopted the face-on thickness of the 
synchrotron emission to be 1 kpc (the typical nonthermal 
disk thickness, e.g. Condon I1992|l . a proton-to-electron 
energy ratio of 100 fPacholczvk ll97()'jl . and a lower cutoff 
energy for relativistic protons of 300 MeV. The effective 
pathlengths were derived from the assumed inclination 
angle. We adopted an inclination of 73° for all regions 
associated with NGC 4038, 55° for NGC 4039 (both from 
LEDA database) and the mean of these values (64°) for 
the whole system and for regions between the galaxies. 

To account for all uncertainties due to a complex three- 
dimensional geometry of the system we repeated the calcu- 
lations for the whole system and the "overlapping region" 
with a mean inclination of 0°. Also, for other regions where 
the value of inclination was uncertain, additional calcula- 
tions were computed for the angle of 0°. For the three 
regions discussed in Sect. 14.1. "21 we performed calculations 
with inclinations determined from our models. All the re- 
sults are presented in Tabled The typical uncertainty in 
the calculated values is about 25-30%, including an un- 
certainty of a factor of two in the nonthermal face-on disk 
thickness, the proton-to-electron energy ratio and in the 
lower proton energy cutoff. 

For the whole system we averaged the radio intensi- 
ties within the area delineated by the intensity level of 
0.3 mJy/b.a. (roughly 1% of the peak intensity) and ob- 
tained a mean total magnetic field strength of i?tot = 
20 ± 5/iG. This is a high value, almost as high as in the 
starburst galaxy M 82 (Klein et al. 1988). The mean reg- 
ular field is Breg = 5 ± 1 /iG and yields a field regularity 
0.25. Compared to other galaxies the mean field regu- 
larity is very low (see Sect. 14 .31) . A smaller face-on non- 
thermal disk thickness would even yield stronger magnetic 
fields. A thickness of 0.5 kpc, as used for NGC 4038/39 by 
Hummel & van der Hulst (^86 ) , and an inclination of 64° 
yield a total field strength of 24 fiG for the whole system. 

The numbers in Tabled show that the total magnetic 
fields in the individual galaxies without the "overlapping 
region" and in the "overlapping region" itself are about 
twice stronger than in normal spirals for which the mean 
value depending on the sample is 8-10 (see Beck et al. 
I1996f) . This field amplification must be a result of interac- 
tion. 

Locally, the strongest magnetic fields of i?tot — 30 ^G 
arc observed outside the galaxy disks, in the southern part 
of the "overlapping region" (C and D in Fig. (Sj), a phe- 
nomenon which was never observed before in any galaxy 
system. This strong magnetic field is likely caused by the 
interaction. In quiescent regions in the northwestern part 
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Table 1. Magnetic fields strength and field regularity in different regions of the merging system of NGC 4038/39. 



Region 


inclination 


Btot 




regularity 




(°) 




f/iG) 




The whole merging system 


64 


20 


5 


0.25 







24 


5 


0.19 


NGC 4038 (without overlapping region) 


55 


21 


5 


0.25 


NGC 4039 (without overlapping region) 


73 


15 


3 


0.19 


Overlapping region as a whole 


64 


27 


4 


0.16 







34 


4 


0.12 


Central peak of NGC 4038 


55 


29 


8 


0.30 


Central peak of NGC 4039 


73 


21 


1 


0.05 


N Dart of the northeastern ridee. modelled in Sect. 14.1.21 


59 


19 


9 


0.52 


(RA= 12''01™53'!0 Dec= -18°51'33", a black circle in Fig.lHJ 










Middle oart of the northeastern ridee. modelled in Sect. 14.1.21 


34 


18 


8 


0.48 


(RA= 12''01™55'!3 Dec= -18°5l'4l", a black star in Fig.lHJ 










SE Dart of the northeastern ridee. modelled in Sect. 14.1.21 


4 


16 


10 


0.88 


(RA= 12''01™57:0 Dec= — 18°52'06", a black cross in Fig. IHJ 










Region A, middle of the northeastern ridge 


55 


22 


13 


0.76 







24 


10 


0.45 


Region B, overlapping region, north 


55 


27 


12 


0.50 







30 


8 


0.29 


Region C, overlapping region, south 


64 


30 


3 


0.12 







37 


2 


0.05 


Region D, overlapping region, south 


64 


32 


7 


0.23 







39 


4 


0.10 


Region E, X-ray ridge 


55 


24 


8 


0.35 







27 


5 


0.20 


Region F, spiral arm in NGC 4038 


55 


23 


9 


0.43 


"Interarm" region in NGC 4038 


55 


22 


5 


0.21 



(RA= 12''01'"5l!4 Dec= -18°52'0", a black triangle in Fig.|Hl 



of NGC 4038 magnetic fields resembles that of a normal 
spiral arm. In region F, which is probably the least dis- 
turbed place of the relic spiral arm located away from 
strong star- forming regions, the magnetic field reaches a 
value of 23 //G (for 1 kpc disk thickness), comparable to 
that observed in the brightest parts of strong spiral arms 
in e.g. NGC 6946 and NGC 1566 (Beck et al. ,1996) . The 
most quiet "interarm" region in NGC 4038 (the last re- 
gion presented in Table ^1 shows i?tot — 22 /xG, much 
higher than for typical interarm regions in normal galax- 
ies. However, this place is not far away from the central 
region of NGC 4038 so that some enhancement of star 
formation is suspected. 

4.3. The regularity and origin of magnetic fields in 
NGC 4038/39 

Particular regions with different field strengths and struc- 
tures can be associated with various phenomena in the 
gas dynamics in NGC 4038/39. Regular magnetic fields 
in NGC 4038/39 are similar to that observed in normal 
spirals and are of order of 5 /iG. However, the field reg- 
ularity is rather low (0.2-0.3 even in the "interarm" re- 
gion). A typical regularity of 0.5 ± 0.1 was found for nor- 
mal galaxies (Buczilowski & Beck lUMl Beck et al.T^M^). 
Our resolution relative to the galaxy size is better than 



in all galaxies used in the quoted papers. In such case 
we should expect even higher degree of field regularity 
(> 0.5) because of lower beam depolarization. This means 
that in the interacting system either strong amplification 
of turbulent fields occurs (see below), or both turbulent 
and regular magnetic fields are generated, but the regular 
fields are much more tangled in the starburst regions than 
in normal spirals. Strong regular magnetic fields (~ 8 fiG) 
are observed close to NGC 4038 center and in its relic 
spiral arm (region F) but still with a low field regularity 
of 0.3-0.4. In the northeastern ridge, far away from star- 
forming regions, the magnetic fields become more ordered 
and reach a strength of > 10 /iC (region A), probably as 
the result of gas shearing motion in the base of the tidal 
tail. 

In the northeastern part of the system (region A) the 
response of magnetic fields to the local gas flows is well 
visible. The magnetic pitch angle change from large values 
(> 40°) in the western and northern part of NGC 4038 to 
small ones in the eastern part, reflecting the alignment of 
regular magnetic field with the tidal tail, probably due to 
the gas stretching motions while a higher flcld regularity 
is suggestive for decaying gas turbulent motions. In the 
northern (region B) and middle part of the "overlapping 
region" magnetic flelds are also aligned with the shape of 
this region possibly adjusting to the gas flows. 
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The enhancement of total magnetic field in the "over- 
lapping region" (C and D) to high values of ~ 30 /xG with 
coherence scale of about 1 kpc (visible in e.g. Fig.|SJ) could 
be produced by the fluctuation dynamo (Subramanian 
I1998|l which may amplify magnetic fields liberated from 
the parent galaxies during the gravitational interaction. 
However, assuming a turbulent magnetic diffusivity of 
about 10^''cm^s~^ (an intermediate value between a typ- 
ical value in the disk and in the corona in galaxies) the 
estimated magnetic diffusion time is about 3 x 10® yr, very 
close to the dynamic scale time of this region of about 
2 X 10® yr (e.g. Mihos et al. 11993(1 . Thus, no amplification 
mechanism by any dynamo action is necessary to sustain 
the magnetic fields in this region if strong enough mag- 
netic fields were already present in the disks. 

Instead we may assume isotropic compressional en- 
hancement of magnetic fields in region C and D accord- 
ing to the approximate scaling with confining gas density 
as B oc p^^^, with 10 of the magnetic field strength 
before compression and the observed value of 30 af- 
ter compression. Current gas density estimated from the 
total mass of gas approximated by the sum of molecu- 
lar gas (from CO observations, Zhu et al. I2003|l and neu- 
tral Hi gas (from Hibbard et al. I2000|l yields a value of 
10.7 cm~^. Using the scaling relation we get a gas density 
of 2.1 cm~'^ before compression. This value is similar to gas 
densities in normal galaxies, e.g. in M31 where gas den- 
sities span a range of 1-6 cm~'^ (Berkhuijsen et al. 11993(1 . 
Hence, compressional enhancement in the "overlapping re- 
gion" of magnetic fields pulled out from the parent disks is 
highly probable. Nevertheless some kind of fluctuation dy- 
namo cannot be excluded, but this needs elaborate MHD 
modelling. 

4.4. Evolution of gas and magnetic field in 
NGC 4038/39 

The coupling of magnetic field with different gas phases 
shows more variety in NGC4038/39 than in normal galax- 
ies where density waves play an ordering role. The pro- 
files of various species along three different directions (dis- 
cussed in Sect. l4.l|l revealed different connections of mag- 
netic fields with cold, warm and hot gas depending on the 
particular location and dominance of various physical pro- 
cesses. The intensity ratios between infrared, X-ray and 
radio emissions yielded additionally a crude estimate of 
the evolutionary stages of the selected places. 

From all the above considerations we propose the fol- 
lowing tentative scenario of the gas and magnetic field 
evolution in NGC 4038/39: 

1. Before the present encounter the northern galaxy was 
gas rich with a strong spiral magnetic field. 

2. The current collision has pulled out the stellar tidal 
tails (in agreement with simulations by Englmaier, 
priv. comm., see also Englmaier et al. 12001(1 . The NW 
part of "inverted 9" structure is possibly a relic of an 
unfolded spiral arm and constitutes the base of an ex- 



panding gas-rich tail extending further to the south. 
The regular magnetic field remains aligned with the 
outer parts of this structure (region A). The magnetic 
field gets more ordered as the turbulence decays in the 
quiet, SE parts of this feature. 

3. Another spiral arm of the original northern galaxy has 
possibly collided with that of the southern galaxy, giv- 
ing rise to the dark cloud complex with a regular field 
aligned along the collision front (regions B, C and D). 
The magnetic field may be already strong in the north- 
ern, cool part of the complex (region B) but there 
is not yet strong star formation able to produce CR 
electrons. Some 5-10 Myr ago a strong star formation 
burst occurred along the NE boundary of the cloud 
(region E), efficiently destroying the regular magnetic 
field and filling this region with a pool of diffuse, hot 
gas. 

4. A young star formation burst has occurred in the 
southern part of the cloud (region D) , being still highly 
obscured in the optical domain, but visible in radio 
and infrared emission. It already managed to disrupt 
the regular magnetic field and to build up exception- 
ally strong total fields by compressional forces. At the 
southern edge of this clouds complex (region C), star 
formation is so young that radio intensity is still weak. 

5. Summary and conclusions 

We present the first study of radio polarized and total 
emission of a pair of spiral galaxies in the process of col- 
lision. The galaxies - NGC 4038/39 - were observed with 
the VLA in the C and D arrays at 8.44 GHz, 4.86 GHz 
and 1.49 GHz. 

We show that the galaxies possess bright, substan- 
tially polarized radio disks with no dominant central 
sources. The thermal fraction is found to be about 50% 
at 10.45 GHz, higher than in normal spirals. The de- 
rived magnetic fields in both galaxies are about two times 
stronger (~ 20 /zG) than in normal spirals probably as the 
result of interaction-enhanced star-formation and stronger 
turbulent motions. A similarly strong mean magnetic field 
is also found along the interface between the two disks 
known as the place of merger-driven extensive starburst. 
In comparison with other galaxies magnetic fields in the 
Antennae have lower regular components relative to the 
total one which is likely due to more strongly tangled fields 
in the vivid star-forming regions. 

Our radio data are compared with existing data in 
H I, Ha, X-rays and in the infrared to study the physical 
conditions in various places of the merger and local in- 
terrelations of magnetic fields with different ISM phases. 
We distinguish several regions in the interacting system 
(depicted as A-F in Fig. IS)) with the following results: 

A. The eastern highly polarized ridge where the gas-rich 
southern tidal tail starts. This region has no counter- 
part in any of star formation-dependent quantities like 
total radio, thermal, nonthermal or infrared intensity. 
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It coincides with the Hi ridge which extends far out 
along the southern tidal tail. The radio spectrum is 
much steeper in this region and still steepens up to 
a > 1.2 towards the tail, which indicates the CR elec- 
tron population aging as they move away from their 
sources in star-forming regions. The magnetic field is 
highly regular (Breg — 10 nG) and stretched along the 
ridge, becoming even more ordered towards the tidal 
tail as turbulent motions decay, it also bends to be- 
come almost parallel to the sky plane. 

B. The cold northern part of a dark cloud complex, consti- 
tuting the "overlapping region" . The radio total, ther- 
mal and nonthermal intensity is rather weak, indicat- 
ing no hidden star formation. We note also that the 
whole region of the dark cloud complex shows a co- 
herent magnetic field structure, probably tracing the 
line of collision between the arms of discussed merging 
spirals. 

C. The hidden star- forming clump at the southern tip of 
the dark cloud complex. The clump is very strong in 
the infrared but deficient in the total radio intensity 
relative to the infrared one. Either the star formation 
burst occurs in rather weak magnetic field around or 
it is so young (< 10^ yr) that not enough CR electrons 
were produced. 

D. The region of hidden star formation in the southern 
part of the dark cloud complex bright in total radio, 
thermal, nonthermal and infrared intensity. However, 
the Ha peak is rather weak due to heavy extinction. 
The region is weakly polarized and most of the mag- 
netic field is random reaching large total values of 
Btot — 30 /iG. This is probably the result of compres- 
sion of the original fields pulled out with gas from par- 
ent galaxies. X-ray emission is relatively weak, possibly 
the starburst is still too young (i.e. < lO'^ yr) to fill this 
volume with hot gas. 

E. The X-ray ridge tracing an older (> 5 x 10^ yr) star 
formation burst leaving a pool of the hot gas and in- 
termediate age star clusters. It has no counterpart in 
total radio intensity (the CR electron population faded 
away) and the regular magnetic field has been effi- 
ciently destroyed by star-forming activity. 

F. The relic of the magnetic spiral of the gas-rich north- 
ern galaxy. The Faraday rotation measure distribution 
shows that this region has a complex three-dimensional 
structure. 

Locally, the strongest magnetic fields of ~ 30 fiG are 
present outside the galaxy disks, in the "overlapping re- 
gion" , what was never observed before in any galaxy sys- 
tem. We also note the large values of the magnetic pitch 
angle in the relic spiral arm of NGC 4038 and the unusual 
"interarm" region with disrupted regular fields. Thus the 
magnetic fields in the Antennae are in structure and prop- 
erties very distinct from what is observed in normal spi- 
rals. 

The presented regions (A-F) reveal different physical 
conditions and evolutionary stages and magnetic fields as- 



sociations with cold, warm and hot gas, depending on 
the particular place and dominance of various physical 
processes. This makes the Antennae a key target for 
studying galactic gas dynamics, the properties of vari- 
ous ISM phases, star- forming processes and their inter- 
relations with the magnetic field. 

We are, however, still far from fully understanding 
the dynamics of NGC 4038/39. Further works involving 
the numerical simulations including the magnetic field 
and high resolution observations of gas kinematics of the 
Antennae in cold, warm neutral, and in ionized gas are 
highly desirable. Especially interesting would be a model 
of gas and magnetic field evolution in the "overlapping 
region" far from the galaxy disks. 
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